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We study the paramagnetic Anderson-Hubbard model using an extension of dynamical mean-field 
theory (DMFT), known as statistical DMFT, that allows us to treat disorder and strong electronic 
correlations on equal footing. An approximate nonlocal Green's function is found for individual 
disorder realizations and then configuration-averaged. We apply this method to two-dimensional 
lattices with up to 1000 sites in the strong disorder limit, where an atomic-limit approximation is 
made for the self-energy. We investigate the scaling of the inverse participation ratio at quarter- 
and half-filling and find a nonmonotonic dependence of the localization length on the interaction 
strength. For strong disorder, we do not find evidence for an insulator-metal transition, and the 
disorder potential becomes unscreened near the Mott transition. Furthermore, strong correlations 
suppress the Altshuler-Aronov density of states anomaly near half-filling. 



I. INTRODUCTION 

The physical properties of interacting disordered mate- 
rials are often qualitatively different from their noninter- 
acting counterparts. For example, it has long been known 
that noninteracting quasiparticles in two-dimensional 
(2D) materials are localized by arbitrarily weak disorder; 
however, there is evidence that interactions can drive an 
insulator-metal transition in 2D.- Similarly, there is a 
growing awareness that disorder can fundamentally al- 
ter the physical properties of interacting systems. This 
arises in a number of transition-metal oxideSf2i'^>^'& where 
the predominantly d-orbital character of the conduction 
electrons results in a large intra-orbital (on-site) Coulomb 
interaction relative to the bandwidth. These materials 
are of interest because they have an interaction-driven 
insulating (Mott-insulating) phase and because of the va- 
riety of exotic phases, such as high temperature super- 
conductivity, which appear near half-filling. However, 
transition metal oxides are typically doped by chemical 
substitution and, with few exceptions, are intrinsically 
disordered. At present, there is little consensus on the 
effects of this disorder, particularly near the transition 
to the Mott-insulating phase. 

Here, we discuss the effects of strong electronic corre- 
lations on disordered two-dimensional (2D) materials via 
a numerical study of the Anderson-Hubbard model, 

H = -t^^cl^Cj^ +^ (f7n,|n,| -|- e.n^) , (1) 

where refers to nearest neighbor lattice sites i and 
J, (T ^tii is the spin index, and fii — -I- fiii where 
fiia = cl^Cia is the local charge density operator. The 
model has four parameters: the kinetic energy t, the 
intra-orbital Coulomb interaction U, the width W of the 
disorder-potential distribution, and the chemical poten- 
tial fi. Disorder is introduced through randomly chosen 
site energies e^, which in this work are box-distributed 
according to P(e,) = W-^e{W/2 - |e,|). 



For U = 0, it is well understood that the single-particle 
eigenstates of Eq. ([1]) are Anderson localized for W > Wc, 
where Wc is the critical disorder and Wc = in two 
and fewer dimensions. For W — 0, and at half-filling 
(ie. n — 1, where n is the charge density), there is a 
critical interaction strength Uc such that the model is 
a gapped Mott insulator for U > Uc and (neglecting 
possible broken-symmetry phases) a strongly-correlated 
metal for U < Uc or for n 7^ 1. There is evidence 
that the Mott transition is fundamentally different in 
the presence of disorder. For example, some work has 
shown that Uc — in clean low-dimensional systems with 
nested Fermi surfaces r'liiSii^ while in the disordered case 
Uc is not only nonzero, but rapidly becomes of order the 
bandwidth as a function of W^^^^^ In fact, it has be- 
come clear through numerous studies that the general U - 
W-n phase diagram is complicated and also potentially 
contains superconducting, antiferromagnetic, and spin 
glass phases.— li^i^ii^ Much of the recent progress has 
been for infinite-dimensional system a"'^'^i"'^^'"'^^i"'^''''"'^^i"'^^'^°'^"'^ 
and the applicability of this work to two and three di- 
mensions is not well established. A recent focus has 
been the extent to which interactions screen the disor- 
der potentia l^^i^^'^^1^^ and whether screening may lead 
to an insulator-metal transition] 12,22,23,24 jj-^ particular, 
some calculations show perfect screening near the Mott 
transitionj ^^'^^ 

Part of the confusion surrounding the Anderson- 
Hubbard phase diagram stems from the variety of 
theoretical approaches that have been applied. Self- 
consistent Hartree-Fock (HF) calculation o^^i^^'^'^ treat 
the disorder potential exactly but do not capture the 
strong-correlation physics of the Mott transition. Dy- 
namical mean-field theory (DMFT) approaches con- 
tain the necessary strong-correlation physics but are 
based on a local approximation that generally pre- 
cludes exact treatment of disorder. A variety of 
coherent-potential-approximation (CPA) and CPA-like 
approximations have been employed in conjunction with 
DMFTpiiii^ii^iiiii^iiai^o While the CPA reproduces some 
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disorder-averaged quantities accurately, eg. the density 
of states (DOS) in the noninteracting hmit, it fails to re- 
produce quantities that depend on explicit knowledge of 
spatial correlations between lattice sites. As a notable 
example, the CPA fails to predict the Altshuler-Aronov 
DOS anomaly which appears at the Fermi energy in dis- 
ordered metals;^ Quantum Monte Carlo^^jSa g^^^^ exact 
diagonalization^"' (ED) methods treat both disorder and 
interactions exactly, but suffer from severe finite-size lim- 
itations, typically generate only equal-time correlations, 
and (in the case of quantum Monte Carlo) suffer from 
the fermion sign problem. 

In this work, we use an extension of DMFT, known as 
statistical DMFTf2£ that incorporates both strong cor- 
relations and an exact treatment of the disorder poten- 
tial. By varying U for a fixed disorder strength, we are 
able to move smoothly from the well-understood weakly- 
correlated regime into the unknown territory of strongly- 
correlated disordered systems. Because our approach re- 
tains spatial correlations between the local self-energy 
at different lattice sites, but can also be applied to rea- 
sonably large lattices in finite dimensions, it provides a 
bridge between the various methods described above. Up 
to now, statistical DMFT has been applied only on a 
Bethe lattice. Here, we work with a two-dimensional 
square lattice and present results for the density of states, 
scaling of the inverse participation ratio, and the screened 
potential. 



II. METHOD 

We focus on paramagnetic solutions on the 2D square 
lattice, for which the noninteracting bandwidth is D = 
8t. On an iV-site lattice, the single-particle Green's func- 
tion can be expressed as an iV x matrix in the site- 
index: 



g,{ij) = [gO(c^)-i - EP(a;)] ' where 



G{uj) = [ojI - t - e - ■E{uj)]- 



(2) 



with I the identity matrix, t the matrix of hopping ampli- 
tudes, e the diagonal matrix of site energies and S(cli) 
the matrix of self-energies. The matrix t has nonzero 



matrix elements ti 



-t for i and j corresponding to 



nearest- neighbor sites. We assume that S(w) is local, 
having only diagonal matrix elements Y,i{Lo). 

The iteration cycle begins with the calculation of G(aj) 
from Eq. 0. For each site z, one defines a Weiss mean 
field a°(w) = [G,j(w)-i -(-E,(w)]-i where Gy (w) are the 
matrix elements of G{uj). As in the conventional DMFT, 
one then solves for the full Green's function Qi{uj) of an 
Anderson impurity whose noninteracting Green's func- 
tion is Q^iuj). The self-energy is then updated according 
to Ef™(w) = e,"(w)"i - g-{Lo)-^ and the iteration cy- 
cle is restarted. In the disorder-free case, this algorithm 
reduces to the conventional DMFT. 

We have used the Hubbard-I (HI) approximation as 
a solver for ^/i(w), for which, in the paramagnetic case. 
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and rii = (hi) is self-consistently determined for each 
site. The HI approximation is the simplest improvement 
over Hartree-Fock (HF) that generates both upper and 
lower Hubbard bands. When either the HF or HI ap- 
proximations are used as a solver, the DMFT result is in 
fact equivalent to the result obtained directly from the 
corresponding approximation. However, the statistical- 
DMFT procedure outlined above, applied to a physical 
lattice, and hence implemented through a matrix inver- 
sion [Eq. ([2|)], is very computationally intensive. In par- 
ticular, it can be difficult to achieve a converged self- 
consistent solution to the N coupled equations for Si(u;). 
This work, aside from being a simple improvement be- 
yond HF, represents an important proof of principle re- 
garding the feasibility of this statistical DMFT approach. 

The fact that we are studying disordered systems bears 
on two important and related issues: (i) the accuracy 
of the HI approximation and (ii) the validity of single- 
site (as opposed to cluster or cellular) DMFT in a finite 
dimensional system. The strengths and weaknesses of 
the HI approximation are well documented in the clean 
limit: it is exact in the atomic limit {U/t ^ 1) but is 
nonconserving and fails to satisfy Luttinger's theoremj^ 
However, it is uniquely effective in low-dimensional dis- 
ordered systems because W/t » 1 corresponds to the 
atomic limit for arbitrary U in systems with finite coor- 
dination number. This is in apparent contradiction to the 
general result that DMFT is only exact in infinite dimen- 
sions, while nonlocal terms in the self-energy, neglected 
in DMFT, can have a dramatic effect on the Mott transi- 
tion in 2Di^i^ In the disordered case, however, we have 
compared our results with ED studies on small clusters^ 
and have found qualitative agreement for W = 12t (ie. 
W — 1. 51?), which is the focus of the current work. By 
pushing the system towards the atomic limit, strong dis- 
order both enhances the accuracy of the HI approxima- 
tion and also reduces the importance of nonlocal terms 
in the self-energy. 



III. RESULTS 

The local density of states (LDOS) is extracted from 
the Green's function as p(ri,w) = —Tr^^lmGa^uj) and 
the density of states (DOS) is p{uj) — N~^^-p{Yi,Lo). 
Figures[TJa)-(c) show the evolution of the DOS as a func- 
tion of U for fixed W at half-filling. There is a transition 
from a single band to a gapped Mott-insulating state at 
a critical interaction Uc ~ 12. 5t. In contrast, the LDOS 
[Figs.[TJe)-(g)] develops a local Mott gap at much smaller 
values of U . We note that, at most sites, the noninter- 
acting LDOS has a single strong resonance near the bare 
site energy [Fig. [IJe)] indicating proximity to the atomic 
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FIG. 1: (color online) Density of States, (a)-(c) Total DOS for 
W = 12t, n = 1 and different U. Results have been averaged 
over 1000 samples and are for an A'^ = 20 x 20 site square 
lattice. ED results, averaged over 10000 samples, for a 4-site 
lattice (red) are also shown in (c). (d) DOS for n — 0.5. Also 
shown (e)-(g) is the LDOS for three arbitrarily-chosen sites 
(site energies indicated in the figure). Panel (e) is for the 
same parameters as (a), (f) as (b), etc. Energies are in units 
oft. 



limit where HI becomes exact. The coordination number 
of the lattice determines the disorder strength required 
to approach the atomic limit. In lattices with infinite 
coordination number, each site couples to a continuum 
of states and the atomic limit is never reached for finite 
W. Thus, somewhat paradoxically, DMFT with the HI 
solver works best, in disordered systems, for lattices with 
low coordination number. 

It is useful to compare the DOS evolution in Fig. [1] 
with existing published work. First, we note that the 
Mott transition occurs at Uc ^ W in both our DMFT 
and ED calculations. This is consistent with QMC re- 
sults for three dimensionsj^ as well as infinite-dimensional 
DMFT results,!- but Uc is somewhat larger than found in 
one dimension.^!*' Our Uc is also significantly larger than 
in unrestricted HF calculations for three dimensionsjii 
however, in HF calculations, the metal-insulator transi- 
tion is of the Slater-type and could therefore be expected 
to respond differently to disorder. 

Previously published results for the DOS near the Mott 
transition in the large-disorder limit, to our knowledge, 
are for infinite dimensions where CPA-like approxima- 
tions can be made. The main distinction between our re- 
sults and those for infinite-dimensions is the quasiparticle 
resonance due to Kondo screening of the local moments 
that appears at the Fermi energy for U ^ Uc in the latter 



case 
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(The peak at uj — /.i for U — 8t comes from 
the overlap of the lower and upper Hubbard bands.) The 
HI solver cannot give such a peak; however, ED calcu- 
lations for small clusters^ also find no resonance at any 
U. In particular. Fig. [^c) shows ED results for parame- 
ters correspon ding to a maximum quasiparticle resonance 
height in Refs. |15li20l |. The difference between our results 
and those for infinite-dimensions might seem unsurpris- 




FIG. 2: (color online) Local density of states for a single 
disorder realisation a,i lj = jj,, W — 12t, for (a) U = 0, (b) 
U = 8t, (c) U = 12t for an = 32 X 32 site lattice and n = 1. 



ing since there is also no resonance in the clean limit in 
2D.^^ However, the reason for this absence appears to be 
different in the clean and disordered systems. The lack 
of a quasiparticle resonance in clean 2D systems has been 
attributed to nonlocal terms in the self-energy,~ terms 
which are not included in HI. These nonlocal terms de- 
scribe short-ranged antiferromagnetic correlations whose 
effect is to suppress Uc to zero in the clean limit.— i^'^^ In 
our work, the absence of a quasiparticle peak stems from 
the calculations being in the atomic limit. That the same 
physics controls the ED results is supported by the fact 
that Uc is nearly identical in HI and ED calculations. It 
therefore appears as if the limitations of the HI solver do 
not hide key physics in the large disorder limit. 

The DOS at quarter-fiUing, Fig. [Ud), does have a 
peak at the Fermi energy, but the origin of this peak 
is unrelated to strong correlations. Paramagnetic HF 
calculations^ have shown that nonlocal charge-density 
correlations lead to a positive DOS anomaly at the 
Fermi level in disordered metals when the interaction is 
zero-range and repulsive. This Altshuler-Aronov DOS 
anomaly is absent at half- filling in our calculations, in 
contradiction with the HF calculations, indicating that 
strong correlations suppress the peak. We will discuss 
this point below. 

In Fig. [21 p(r, fi) is plotted for different values of U for 
a particular disorder configuration. At t/ = 0, sites which 
have significant spectral weight at lo = fj, are typically iso- 
lated from one another, consistent with electrons being 
Anderson-localized. The LDOS is more homogeneous for 
U ^ 8t than for U = 0, consistent with an interaction- 
driven delocalizing effect. However, at U ^ 12t the 
LDOS is again highly inhomogeneous. 

The inverse participation ratio (IPR) 



l2iL0,N) 



(4) 



provides a quantitative measure of the inhomogeneity of 
p{r, Lo). The IPR can also be used to distinguish extended 
and locahzed states since limjv_oo -f2('^, ^) — for the 
former and is nonzero for the latter. It is important to 
note that the frequency lo used in the calculation of G{lo) 
in Eq. ^ contains, by necessity, a small imaginary com- 
ponent 17. Scaling quantities such as the IPR that are 
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FIG. 3: (color online) IPR scaling with system size for W = 
12t. Shown are (a) the dependence of the IPR on 7 for U = 0, 
(b) the scaling of the IPR with 7 — it/N and n = 1, and 
the dependence of the localization length (defined in text) 
on U a,t (c) n — 1 and (d) n = 0.5. All curves are for uj — 
fj,. Anderson-insulating (AI) and Mott-insulating (MI) phases 
are indicated in (c). 
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FIG. 4: (color online) (a),(b) Local charge density and (c),(d) 
corresponding screened potential as a function of site energy. 
The screened potential in the t = limit (red dashed) is also 
shown in (d). The relative variance of Vi (circles) is shown 
(inset) along with the t = result (dashed curve). 



derived from G(w) generally depend on 7<^ In our scal- 
ing calculations, we have taken 7 oc 1/N such that the 
ratio of 7 to the level spacing remains constant. As shown 
in Fig. [S^a) for the noninteracting case, the effect of 7 is 
to reduce hi^^, N), such that our results represent a lower 
bound on the true (7 — 0) IPR. The IPR scaling at half- 
filling is shown in Fig. ^h) for uj — jj.. For each value 
of U, we extrapolate a limiting value /2(/^, 00). We then 
define a localization length ^ — /2(/i, 00)^^^'^, where d is 
the dimension of the system. For finite 7, I2i^J■, oo)"^/'' 
gives an upper bound for ^. 

FigureslSIc) and (d) illustrate the effect of strong corre- 
lations on localization. While the HI results agree closely 
with self-consistent HF calculations for small U, the dis- 
crepancy between the methods grows as U is increased. 
In the weakly-correlated small-?/ regime, ^ grows with 
U, consistent with increased screening of the impurity 
potential . ^^'^'^ For large U, however, ^ is a decreasing 
function of U and is smaller near the Mott transition 
than in the U = case. As we discuss below, this can 
be partially attributed to a decrease in screening due 
to strong correlations, although screening no longer pro- 
vides a complete framework for understanding the evo- 
lution of ^. 

To understand better the localizing effect of strong cor- 
relations, we define a screened potentia l^^'^°i^^ Vi based 
on the HF site-energy: V^ = + [7^^,^ Plots of both 
Ui and Vi as a function of are shown in Fig. 3) In the 
weakly-correlated limit {U — t), rii is approximately lin- 
ear in ei over a wide range and Vi « ei(l — Uxu) with 
Xii = —dui/dei. Since xu depends only weakly on J7, Vi 
is a decreasing function of U. In the strongly-correlated 
limit (U = 8t), Ui is a nonlinear function of e^. For 



ej| < J7/2, rii « 1 and Xu — ^^ such that these sites are 
unscreened (ie. Vi = ei). There are, therefore, two limits 
in which screening is small at half-filling: U <^W where 
the interaction is too weak to screen the impurity poten- 
tial, and U ^ W where strong correlations enforce single- 
occupancy of most sites. We note that the absence of an 
Altshuler-Aronov DOS anomaly near half-filling can be 
understood in this context: strong correlations suppress 
the local response of the charge density to the impurity 
potential. 

A measure of the screening is given by the relative 
variance AV'^, defined as the variance of Vi divided by 
the variance of e^. Our numerical results. Fig. [4l show 
that AV'^ is a nonmonotonic function of U obtaining a 
minimum at [/ w W/2 and approaching AV'^ = 1 for 
t/ — > and U > W. At a qualitative level, this is consis- 
tent with the nonmonotonic dependence of ^ on [/, since 
one expects ^ to be large when AV'^ is small. There 
are, however, quantitative discrepancies which show that 
AV'^ does not tell the whole story. First, ^ does not ob- 
tain its maximum a,t U — W/2, where AV'^ obtains its 
minimum. Second, ^ is smaller at large U than at U = 0, 
indicating that states near the Mott transition are more 
strongly localized than at U = 0. 

Our results for the IPR are consistent with re- 
cent quantum Monte Carlo calculations of the dc 
conductivit y^^'^^ which find a similar nonmonotonic de- 
pendence on U. The results for the screened poten- 
tial, however, are inconsistent with infinite-dimensional 
DMFT results J^i^o In Ref. impurities are found to 
be perfectly screened (ie. AV'^ 0) at the Mott tran- 
sition, which would correspond to a divergent £ in our 
calculations. Near the Mott transition, Ref. [21| found 
that sites with \ei\ < U/2 are perfectly screened. Both of 
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these results are opposite to what we have found here. 

There are two important distinctions between our cal- 
culation and those of Refs. T^Hlj, both of which con- 
tribute to these opposing results on screening. First, 
Refs. [illll use an effective medium approach for the 
disorder potential that results in metallic behavior for 
small U. This is reasonable in high dimensions where An- 
derson localization only occurs for strong disorder. Our 
exact treatment of the disorder potential, on the other 
hand, allows us to describe the Anderson localized phase 
that occurs in 2D for small U. Whereas the LDOS is con- 
tinuous and relatively uniform in the metallic case, it is 
inhomogeneous and dominated by small numbers of res- 
onances in the Anderson insulating phase [cf. Fig. [TJf)]. 
The local charge susceptibility xa is suppressed for sites 
with small LDOS at the Fermi level, and screening in 
the Anderson-insulating phase is consequently expected 
to be less than in the metallic phase. The second dis- 
tinction is the quasiparticle resonance which arises in the 
infinite-dimensional case, but does not occur in our cal- 
culations. This resonance is a key factor in the perfect 
screening found in Refs. [19II2H . 



In conclusion, we have studied the 2D Anderson- 
Hubbard model at half- and quarter-filling, in the limit 
of large disorder using statistical DMFT. We have cal- 
culated the localization length ^ from the inverse par- 
ticipation ratio, and find that it varies nonmonotoni- 
cally with the strength of the interaction: at small U, 
the interaction screens the impurity potential, but at 
large U strong correlations reduce the screening. As a 
consequence, the Altshuler-Aronov DOS anomaly is sup- 
pressed at half-filling. For strong disorder, we find no ev- 
idence for an insulator-metal transition nor for enhanced 
screening near the Mott transition. 
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